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We propose to use the coincidence method of Ma to measure an entropy of the system created in heavy ion 
collisions. Moreover we estimate, in a simple model, the values of parameters for which the thermodynamical 
behaviour sets in. 
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Existence of thermodynamical equilibrium in 
heavy ion collisions is an important question. 
Many phenomenological models critically depend 
on this assumption, and consequently there ex- 
ists an ongoing theoretical debate of this problem. 
Alternatively one should develop methods which 
would allow to answer this question directly by 
the exeriments. For example how to check the 
saddle point relation 
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which is central in a thermodynamical descrip- 
tion. Out of four observables entering above 
equation (energy E, multiplicity n and tempera- 
ture T), the entropy S is most difficult to directly 
measure in experiment. In this talk the recent 
proposal to determine entropy experimentally is 
discussed [Q. 

The Boltzmann relation S(E,n) = logT(E,n) 
reduces the problem to measuring the density of 
states T(E, n) of a system with given energy and 
multiplicity. To this end we have proposed to em- 
ploy the coincidence method advocated by Ma Q 
[]. Suppose that the phase-space is divided into 
multidimensional cells. Suppose furthermore that 
our system occupies T cells (with a uniform prob- 
ability). Each cell represents a different state of 
the system. Our problem is to calculate T. Let 
us select randomly Af configurations of the sys- 
tem. The main advantage of this approach is that 
Af ~ VT <C r is sufficient. These configurations 
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2 A simple variant of this technique is used in common life 
by fishermen to estimate the nubmer of fish in a lake. We 
thank K. Zalewski for this comment. 



occupy some cells. The average occupation num- 
ber of a cell is Af/T <C 1. Under this condition, 
the average number of pairs in the same cell is 
(Af/T) 2 The total number of coincidences N c is 
the sum over all cells N c = T(Af/T) 2 , hence 



r = 



Nt 



(2) 



where Nt ~ Af 2 is the total number of pairs. 
Generalization for nonequivalent configurations is 
also available Jl]||. 

To define meaningfully a coincidence of two 
states, we discretize particle momenta pi = a rii 
with some discretization scale a. As usual a 
should be chosen such that computations are 
practically feasible and at the same time repro- 
duce continuum physics. 

We have tested the method in the case of the 
classical gas of noninteracing, nonrelativistic par- 
ticles with mass m in d dimensions. The dis- 
cretized expression for the number of momentum 
states of N degrees of freedom with the total en- 
ergy E = Ma 2 /2m reads 



r(M, N) 
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n 1 ,...n N ,nf + ...nj 



= M 



where the momenta pi = an, as above. 

Numbers T(M, N) satisfy simple recursion rela- 
tion which will be used to calculate them exactly. 
These will serve to benchmark the performance 
of our Monte Carlo. 

Two questions were addressed: a) how large 
numbers T can be reproduced by the coincidence 
method with present computers and b) how well 
is the thermodynamic limit, M, N — * oo, M/N = 
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e-fixed, 

i lo g r(M, N) - -[log (e) + log (2tt) + 1]. (4) 

approximated within available window of M and 
N. 

A sample of runs is summarized in Table 1 . In- 
stead of generating Af configurations of N integer- 
valued momenta {ni, n,2, ■ ■ ■ , njvjfc, k = 1, ...,Af. 
we have uniquelly labelled each, fc-th say, config- 
uration by an integer index Ik, k — 1, .., T(M, N). 
Consequently each Monte Carlo run consisted of 
a generation of a sample of Af integer indices, 
{h,I 2 ,...,Ix)A < h < T(M,N),k= l,...,Af, 
uniformly distributed in the whole space of avail- 
able states. Then we counted the number of coin- 
cidences N c , i.e., the number of pairs (Ij, Ik) such 
that Ij = Ik- The estimate for the number of all 
states (column 5 of Table 1) is then 

f=Af(Af-l)/N c . (5) 

Multinomial nature of the above process allows 
for simple calculation of the distribution of N c In 
particular, the dispersion of N c reads (Af >> 1) 

CT 2 [A C ]=2< A c >(l + 2^ + ^), (6) 

with < N c >= A" 2 /T, cf.(g). This gives for the 
relative error of the determination of T after Af 
trials y / CT 5 [fJ/r = V2T/Af and for the estimate 

of the error a[T]/f = J2/N c . The last estimate 
is quoted in column 6-th while the actual rela- 
tive deviation is in the last column. It is evident 
from these formulas that the coincidence method 
works for much smaller number of trials (~ VT) 
than the standard approach which measures av- 
erage occupation of a single state. The estimated 
error is steadily decreasing like 1/Af and actual 
deviation follows the suit albeit with some fluc- 
tuations. It is interesting to note that the errors 
decrease as a number of trials and not as 1/ V77, 
sec 0] for more details. Substantial improve- 
ment in the performance can be achieved if naive 
counting of pairs is replaced by " binwise" count- 
ing, i.e. a set of generated indeces {I\, . . . , Ijj} 
is first sorted. This trick reduces the computing 
effort from 0{N 2 ) to O {Af log AO- We find that 



the Monte Carlo results are well under contol and 
show that the method is quite reliable. With the 
sorting trick it is practical for T of the order of 
10 10 . We will discuss now if this is sufficient to 
see the onset of thermodynamic properties. 

Figure 1 shows the entropy density as a func- 
tion of the scaling variable e = M/N. Statistical 
errors of MC results (and the deviation from the 
exact discrete values given by T(M, N)) are much 
smaller that the size of symbols. Exact values are 
very close to the continuum expressions for the 
volume of the nonrelativistic phase space (solid 
lines) . Considered as a function of e and N they 
obviously show a substantial A-dependence. The 
N varies from 8 (lowest curve) to 24 in this plot. 
On the other hand, the deviation from the ulti- 
mate scaling limit, (Eq.(^), the uppermost curve), 
is around 30% in the worst case (N=8,M=30). 
With N starting from 12, deviations from the in- 
finite system are smaller than 20%. Note that A 
denotes the number of degrees of freedom, which 
in d space dimensions corresponds to N/d parti- 
cles. 
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Figure 1. Entropy density s — logr(M, N) vs. 
e. Black symbols represent our Monte Carlo re- 
sults for N=8 (diamonds), 12 (circles), 16 (boxes) 
and 24 (a triangle). 



Second, we test the saddle point relation ([[]) 
with the temperature eliminated with the aid of 
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800 000 


10 


6.39999io - 


-10 


0.44 


0.24 


31 
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1 600 000 


42 


6.09523io - 


-10 


0.22 


0.18 






3 000 000 


188 


4.78723io - 


-10 


0.10 


0.08 


r 








51 795 303 424 










750 000 


12 


4.68749i - 


-10 


0.40 


0.13 


13 


24 


1 500 000 


74 


3.04054io - 


-10 


0.16 


0.27 






3 000 000 


256 


3.51562io- 


-10 


0.09 


0.15 


r 








41 469 483 552 







Table 1 



Monte Carlo results for f (M, N) (col. 5) . The third and fourth column give the number of generated 
configurations AT, and the number of observed coincidences N c . In the last two columns we quote the 
Monte Carlo estimate of the relative error, see the text, and the actual relative deviation S/T = \T — T\/T 
from the exact value T also quoted in the Table. 



the equipartition relation E/N 
cretization it reads 



log 



T(M + 1,N) 
T(M, N) 



N 
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T/2. Upon dis- 
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This equation is tested in Fig. 2, where a half 
of the inverse of the left hand side, as obtained 
from simulations, is plotted as a function of e. 
Solid line represents the right hand side 0. Simi- 
larly to the previous case agreement is very good 
for N > 12. It was necessary to reduce MC er- 
rors to the level of l%-3% in order to achieve this 
agreement. Of couse this test is much more sensi- 
tive than the previous one since it requires precise 
measurement of the derivatives. 

To conclude, the coincidence method is satis- 
factory in practice for the number of degrees of 
freedom (or number of cells) up to ~ 30. This is 
sufficient to see the signatures of the thermal equi- 
librium. For more than 12 degrees of fredom the 
scaling of the entropy density is confirmed with 
the accuracy better than 20% . The saddle point 
relation, coupled with the equipartition principle, 
dS/dE = l/2e is also very well reproduced. 

This work is supported in part by the Pol- 
ish Committee for Scientific Research under the 
grants no. 2P03B 08614 and 2P03B 04412. 



Figure 2. Testing the relation (0). Symbols as 
in Fig.l, solid line corresponds to the thermody- 
namical lmit. 
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